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In this paper we review some properties for the evolving wormhole solution of Einstein equations 
coupled with nonlinear electrodynamics. We integrate the geodesic equations in the effective geom- 
etry obeyed by photons; we check out the weak field limit and find the traversability conditions. 
Then we analyze the case when the lagrangian depends on two electromagnetic invariants and it 
00 ' turns out that there is not a more general solution within the assumed geometry. 
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DhI I. INTRODUCTION 

<: 

^+ . Recently, the interest in worniholes has increased because of the possibihty of intersteUar travel or future time travel 

CnJ ' to past world. Since the formulation of Einstein's equations, several features of these solutions were addressed, such as 

Einstein- Rosen bridges [T]; Wheeler ^ analyzed this kind of solutions and coined the name of wormholes. The most 

fj ' important and popular contribution to these solutions was given by Morris and Thorne [3|, they were the first who 

Q^l made a complete and detailed investigation in this field. Wormholes are solutions to Einstein's field equations which 

I ■ represent two connected universes or a connection between two distant regions of universe. For realistic models of 

{^Tj] traversable wormholes, one should address additional engineering issues such as tidal effects. People is skeptic about 

I— ', them because they contradict some commonly accepted reasonable energy conditions, for example, their structure 

requires matter with negative energy density as a source, i. e. exotic matter [s!]. A full revision on wormholes can be 

►^ , consulted in [J| and 0. 

^^ ' Among an enormous variety of wormhole solutions, there are some kind of spherically symmetric wormholes that 

^^ J evolve with time. This type of wormholes have been extensively analyzed. Roman ^] explored the possibility that 
0^ inflation might provide a mechanism to enlarge tiny wormholes to macroscopic size and investigated the possibility of 
CO ■ avoiding the violation of the energy conditions in the process. Further evolving dynamic wormhole geometries were 
■«H- ■ analyzed, considering specific cases [7|- Q; and also was considered the evolution of a wormhole model in an FRW 
^3 ■ background [W^ . 

00 , On the other hand, nonlinear electrodynamics has recently been applied in several branches in physics, namely, 

^^ ' as effective theories at different levels of string/M-theory [ll|, cosmological models [12], black holes [IJl- [lj| and 
^ in wormhole physics [l5|, among others. In this context, in a recent paper it was shown that (2 -t- 1) and (3 
+ l)-dimensional static spherically symmetric as well as stationary axisymmetric traversable wormholes cannot be 
supported by nonlinear electrodynamics [16]. However, it was found an evolving wormhole solution within nonlinear 
S I electrodynamics |17| . where the lagrangian depends in nonlinear way on one electromagnetic invariant. 
■ " ' In this paper we explore some properties of the evolving wormhole geometry in the context of nonlinear electrody- 

namics. This paper is outlined in the following manner. In section II, we describe briefly (3-|-l)-dimensional evolving 
spherically symmetric wormholes coupled with nonlinear electrodynamics. Section III is devoted to the study of 
wormhole geometry in the aspects of: traversability conditions, linear limit and photon trajectories. In Section IV 
we analyze the case when the lagrangian depends on the two electromagnetic invariants, attempting to find a more 
general solution. In section V some conclusions are given. We shall use geometrized units, i.e., G = c = 1, throughout 
this work. 

II. (3 + 1) DIMENSIONAL EVOLVING WORMHOLE SOLUTION 
A. Action and spacetime metric 

The action of (3 + l)-dimensional general relativity coupled to nonlinear electrodynamics is given by 
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where R is the Ricci scalar. L(F) is the lagrangian, depending in nonUnear form on a single invariant F given by 
F = F^'^F^y/4, [1^, where F^i, is the electromagnetic tensor. Note that in Einstein-Maxwell theory the lagrangian 
takes the form L{F) = —F/Att, however, we shall consider more general electromagnetic lagrangians. 

Varying the action with respect to the gravitational field provides the Einstein field equations G^,y = SnTf^^, with 
the stress-energy tensor given by 



T^i, — 5^1/ L{F) — F^aFy"" Lp . 



(2) 



where Lp = dL/dF. The variation with respect to the electromagnetic potential A^, where F^^ — A^^, — A^i^^, yields 
the electromagnetic field equations 



[F^"" Lp). 
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The spacetime metric representing a dynamic spherically symmetric (3 -I- f )-dimensional wormhole, which is con- 
formally related to the static wormhole geometry Q, takes the form 



ds^ = Q?-{t) 



52*(-)di2 



dr^ 



1 - h{r)/r 



r^{d0^ 



(4) 



where $ and b are functions of r, and Q — f2(t) is the conformal factor that is finite and positive definite throughout 
the domain of t. $ is the redshift function, and b{r) is the shape function [3|. We shall also assume that these 
functions satisfy all the conditions required for a wormhole solution, namely, <I>(r) is finite everywhere in order to 
avoid the presence of event horizons; b{r)/r < 1, with t'(ro) = tq at the throat; as well as the fulfilment of the flaring 
out condition, b — b'^ > 0, with b'{ro) < 1 at the throat. 

Now, taking into account the metric ^ the electromagnetic tensor compatible with the assumed symmetries is 
given by 



F,, = £;(x") (<5* ^S^ , ~ S^ /* J + B{x") {6' ^5^ , - 5^ ^6' J , 

being the non-zero components the following: Ftr = —Frt = E, the electric field; and Fg^ = —F^ 
field. The invariant F takes the form 



F = 



1 



1 _ _ I e-2* E^ - 



B^ 



(5) 
_B, the magnetic 

(6) 



B. Einstein field equations 

We recall that from [l7[ > the solution requirement $ = must be imposed, so using an orthonormal reference frame 
the non zero components of the Einstein tensor reduce to 
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where overdot denotes a derivative with respect to the time coordinate, t, and the prime a derivative with respect to 
r. Analogously, the nonzero components of the stress energy tensor, from Eq. ^, take the form 



Tu = -L- 






J- nn 



L 






Til 



n^r^ sin^ 



-B^Lf 



(8) 

(9) 

(10) 



It is clear that T^^ = 



'Tff and using Einstein's field equations, the following relation is obtained 
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In Eq. (jlip each side depends only on one variable, therefore a solution is found by separating variables: 



b{r)=T[l-a\r^-Tl)\ , 



(12) 



VL{t) = 
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where a is the separation constant, Ci and C2 are constants of integration. If Ci = C2, ^i is singular at t = 0; 
therefore for a > we need to impose Ci > C2 and if a < it is required that Ci < C2, otherwise the conformal 
factor becomes singular somewhere in the domain of t. The conformal function r2(i) — > as i ^ 00, which reflects a 
contracting wormholc solution. It is important to observe that f2(i) — > 00 as t — > io = Q;^^log(±Cj" VC\C2) showing 
a time singularity that must be avoided. In the spirit of [17| we define the dimensionless parameter Rq — aro, so that 
the shape function is given by 



b{r) =rl\-Rl 



ro 



(14) 



A fundamental condition to be a reliable wormhole solution is imposed, that b{r) > [1^. Thus, the range of r is 
ro < r < a = ro\/l + l/i?§; the latter may be arbitrarily large by taking _Ro ^ 0. If a » tq, i.e., Rq ~ ro/a <ti 1, one 
may have an arbitrarily large wormhole. 



C. Electromagnetic field equations 

Solving the electromagnetic field Eq. Q, together with 



(*^P^) 0, 



(15) 



that can be deduced from Bianchi identities, where the asterisk denotes the Hodge dual [i^| , we obtain the restrictions 
that Ftr — —Ftr ~ E{t, r), Fg^ — —F^q = B{9), Lp = Lp{t, r). Thus the magnetic field is given by 



B{e) ^ q^, sine , 



(16) 



where qm is a constant related to magnetic charge. 
Furthermore, from Eqs. ([8]) and pO)) . we obtain 
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The regular solution appears when we set E — 0. Using Eq. (flT)) we obtain 
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and taking into account Eqs. ([7]) and ^, the lagrangian is given by 



SttH^ 
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These equations, together with B = qy^smO, E = 0, F ~ gj^/(2il"'r"') and functions p2)l and ([13]) give a wormhole 
solution well behaved at the throat, with finite fields. This result is in close relationship to the regular magnetic 
black holes coupled to nonlinear electrodynamics found in [ij]. We remind that for static solutions with nonlinear 
electrodynamics, the null energy condition is not violated, fact that forbids the possibility of wormholes; however for 
the evolving wormhole under study this does not apply and moreover, the weak energy condition is satisfied. 

III. ANALYSIS OF THE E^O SOLUTION 

In what follows we address mainly four aspects of the solution, namely: flaring out condition and embedded diagram, 
traversability conditions, the weak field limit and photon trajectories. 

A. Flaring Out Condition and Embedded diagram 

As can be found in [3|, [j|, an imperative feature to appear in any wormhole solution is the so called flare out 
condition once the embedded diagram is given. The embedding is obtained when one considers the wormhole metric 
for an equatorial slice 9 = 7r/2 and in a fixed time t, embedding then in a flat three dimensional Euclidean space, 
ds^ = dz'^ + dr'^ + r^d(f>'^. 

In our case the flare out condition requires ([i3|) that 



at or near the throat. Using Eqs. JT 



>0 



d'^fj-z) _ 1 d'^r (z) _ 1 b-b' 
dz2 - n (i) dz2 ~ n (t) 262 

and P^ the flare out condition can be casted as 
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If a > 0, then we need to impose Cie"* > C2e~"*, and if a < and Cie"* < C2e^"*, otherwise the conformal 
factor becomes negative somewhere along the domain of t, violating the condition Q (t) > 0. So in general the flare 
out condition is fulfilled by the solutions p4| and (fT3l) . 

Another feature that must be explored for a wormhole solution is the one related to the form of the function z (r) 
from the embedded diagram [3| , [J| , this function can be obtained integrating the relationship 



dz (r) dz (r) 



± 



d/r dr 

This equation can be rewritten using equation (|14p to obtain 
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from which we would be able to obtain the solution for z (r) as 



zir) = ±n{t)z{r) ^±n{t) 
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In Fig. [U [20| it is shown the embedded wormhole z{r) for an arbitrary constant time. We can observe that the 
behavior of the derivative is highly positive near the throat ro, as was expected for a wormhole solution, and as r 
grows far from the throat, the value of the derivative ^ becomes smaller but positive, indicating that the tangent 
lines for z(r) become nearly horizontal. This behavior grant us a shape for the throat very similar to the typical 
behavior expected for a wormhole solution [3] , [J| . 
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FIG. 1: The embedded wormhole z{r) for an arbitrary constant time is shown for distinct values of the parameter a. 
continuous curve corresponds to a = 0.2 while the dotted one is for a = 0.3 



The 



B. Traversability Conditions 



We explore now the possibility that a human traveler traverses an evolving wormhole like the one presented here. 
Leaving aside the study of the stability of the solution, we delve into the analysis of the tidal gravitational forces that 
an infalling radial observer must bear while traversing. We take as our traversability criteria the magnitude of the 
tidal forces that an observer can stand during the trip: it must not exceed the forces due to Earth gravity, much on 
the spirit of 21| . To simplify our calculations we shall work in a proper reference frame given by 



-6 = 764 T lf3ef , ej = Tief + 7/36^ , 62 = 6^, eg = e : , 



(25) 



where 7 = , and f3 — v/c, v is the velocity of the radial motion 9]. Then the traversability criteria can be 

stated as 
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where g^ = 9.81m/s^ is Earth gravity; these conditions can be fulfilled with the appropriate selection of the values 
for Ci and C2, rendering this solution in a traversable one. We point out the fact that Eq. (P7)) does not depend on 
the traverse velocity v, however it is not in contradiction with the results given in [J|, Chap. 13.1.1., related to the 
dependence of tidal forces on the traveler velocity, due to the fact that the condition Rgf§f — —Rg^gf. does not lead to 
the presence of a horizon, instead, it take us to equation (jlip from which the solution is obtained. 

Since in the studied case the redshift has been set equal to zero, the radial tides are zero, then if some danger exists 
for the traveler it is in the nonzero accelerations while traversing the wormhole. In what follows we show that there 
is not tidally induced shear. 

The Qij matrix that describes tidal forces is defined by 

Qij = -Ri.a.pivirv'ivjrv^, (28) 

where rjj, J = 1, 2, 3 is an orthonormal triad orthogonal to the four-velocity V". Let us consider the traveler's motion 
is non radial with a velocity vector given by 

V^" ^{j,j(3cosip,0,-ff3sm^P). (29) 

In the case we are analyzing, the tidal forces are described with the entrances 

Qii = -i?friv cos^ V - Ri^i^ sin^ ^ = C1C2 , (30) 

O22 = -7\Rifif sin' V' + Rm<}, cos' ^ + P'Rf4,f<p) = C1C2 , (31) 

Q33 = -^'iRi^i^+P'R,^,^cos'^ + f5'Rg^g^sm'^) (32) 



and 



-l'hC,C2 + C,C2p' cos2 V + ^^^^[(Cie"* + C2e-"*)2] 
+ 4^(1 - «'(^' - r^mCie'^' - C2e-'^r]} , 



Qi2 = -liRifif - Ri^i$) sin V' cos V' = , (33) 



h3 = Q23 = 0. (34) 



Therefore, the tidal forces do not diverge, furthermore, since Q12 — 0, there is no presence of tidally induced shear, 
this means that the traveler will pass safe through the wormhole as far as tidal forces concern. 

C. Maxwellian Limit. 

It is important that any nonlinear theory recovers the form and results of the corresponding linear theory, nonlinear 
electrodynamics is not an exception |18j . However there are cases where the Maxwellian limit is not recovered J22| . 

To carry on the analysis of the Maxwell limit, it is useful to rewrite the main functions F, L and Lp in terms of 
the coordinates r and i, 

L - -^ [1 - "'(3-' - -o)][^^i^ ^^^]' + 3a2(Cie- + C^e'-)^ , (36) 

327ra^ r I 



The result of the analysis of the solution is the absence of a complete Maxwellian limit, i.e., when r ^ oo we get 
that F ^- 0, L ^ ~'^^%r^ ^'^'^ ^F ~^ °°- ^^^ ^^'^^ ^^^^ ^^"^ L lagrangian L does not go as F, and instead goes to a 
constant value shows that in this limit the spacetime is not asymptotically flat due to a constant remanent energy. 
Another impressive result is the corresponding to Lp which goes to infinity, showing a similar behavior to an ideal 
conductor ^23)] . 

Alternatively, as the solution is only valid for the region vq < r < a = ro\/l + 1//?^, and the Maxwellian limit 
would be attained as r — > oo, the requirement of the presence of a region with Maxwellian limit losses strength. 
Furthermore, it does make sense the absence of such a limit since for weak fields the energy conditions are always 
fulfilled, a fact that is contrary to the existence of wormholes. 

D. Light rays in the NLED effective geometry 

In nonlinear electrodynamics photons do not propagate along null geodesies of the background geometr y, in stead 
they propagate along null geodesies of an effective geometry which depends on the nonlinear theory [18| , [2J| . The 
discontinuities of the field propagate obeying the equation for the characteristic surfaces S. For a curved spacetime 
the equation for the characteristic surfaces is 

g^'^S^^S^, = (38) 

And when nonlinear electrodynamics is involved, the corresponding "eikonal" equation for the propagation vectors 
kt" is 

{LFg''" - ALffF'"'F^) k^^k^, = g^^fk^^k^, - (39) 

In the orthonormal tetrad, the expression of g^ff in terms of the electromagnetic field tensor T^'^ also shows clearly 
that the modification in the trajectories is due to nonlinear electromagnetic field: 

''ff - (^^ + ^) ''"' + T7^"' ('') 

given in orthonormal coordinates fii), with 77'"^ — diag[— 1, 1, 1, 1]. For the case of the evolving wormholc with nonlinear 
electrodynamics (JU with P^ and (fT5|) . the effective geometry is given by 

4/ - -f^i^ + si^) =-<?:}/, (41) 



gl^f = Lf + 2L^ 2FLff = g^J^, (42) 

where we have taken into account the energy-momentum tensor given by 

diag(T*^T^^^^^^'^'^) 

with 

T" = -L = -^''^ (43) 

T^^ = L- 2FLf = T^^, (44) 

that are Eqs. (18)-(20) in fV^. Moreover, the expression for the invariant F, if £^ = is 
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In the effective metric the geodesic equations are 
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(45) 
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(47) 
(48) 

(49) 
(50) 
(51) 
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where r is the affine parameter that generates the geodesic trajectories r{T), t{T), 9{t), ^{t). A is a conserved quantity 
related to the existence of Killing vectors. The Eq. ([5T|) is derived from the line element and S — 1. 0, —1 for timelike, 
null and spacelike geodesies, respectively. 

Substituting ((^)2 - (^)2) from Eq.l^ in Eqs. ^ and ([50]), for 5 = it is obtained 
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The constant of motion A makes easier the integration of the geodesic equations. Considering the expressions for 
Lf, L,F in terms of r and t, straightforward algebra leads to the expressions 
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Eqs. ((55)) can be integrated by multiplying the first equation by 2{dt)/{dT) and the second by 2{dr)/{dT), so they 
acquire the form of 
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that allows one to obtain the trajectories r(i) of photons as 
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where Ki and K2 are constants. Substituting f(r,t) and g{r,t) from (|56|) and simplifying we get 



(59) 
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(60) 



where B = C1/C2, -B < 1 for the square root be well defined; Ci and C2 can be adjusted to fulfil the traversability 
conditions as well. 

The graphics [20| of the trajectories r(i) are shown in Fig[2] for different values of a, with Ki = K2 — 0, showing 
the influence of the magnetic field variation. FigO shows r{t) when Ki and K2 are not zero. It can be observed in 
both plots the nonvanishing throat, as well as the non-flatness of the spacetime (associated to the non-Maxwellian 
limit). The effective spacetime as "seen" by photons does not carry as much restrictions as the mere existence of the 
wormhole. 




0.7 2.1 



3.5 4.9 



FIG. 2: Null geodesies rit) in the effective geometry of the nonlinear electromagnetic field. In these plots Ki = K2 = and the 
static case is for a = 0. The non-flatness of the spacetime is apparent, as the geodesies do not tend to infinity for large times 



IV. EVOLVING WORMHOLE SOLUTION WITH TWO INVARIANTS F AND G 



The previous analyzed wormhole corresponds to a lagrangian that depends only on one of the two electromagnetic 
invariants, F = F'^'^F^^/A. On the search for a more general solution we explore in this section the case of a lagrangian 
depending in nonlinear way on both invariants. 
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FIG. 3: Trajectories of photons r{t) in the effective geometry. The continuous curve corresponds to Ki = 1.8, K2 ~ 0.5, while 
the dotted one is for Ki = 0.5, K2 = —0.66, in both a — 0.8. The nonvanishing of the constants of movement enhances the 
action of the magnetic field and smooths the non-flatness of the spacetime 

Let us consider the action of (3 + 1)— dimensional general relativity coupled to nonlinear electrodynamics given by 



5 = 



~9[^ + L{F,G)]d% 



(61) 



where R is the Ricci scalar. L{F, G) is a gauge-invariant electromagnetic lagrangian depending on the invariants 
F = \F^'^F^y and G = jF'^'^F^^ where F^^ is the electromagnetic tensor and F^^ = ^e^^a/3-F"'^ its dual. 

Varying the action with respect to the gravitational field provides the Einstein field equations G^jy ~ SttT^^, with 
the stress-energy tensor given by 



T^i^ — g^_LvL — F^sF^, Lp — F^sF^ Lq, 



(62) 



^^ and Lg = j^ are functions of r and t only, in the same form L — L{t, r). The variation with respect 



where Lp - gp ^..^ ^^ - m 

to the electromagnetic potential A^, yields the electromagnetic field equations 



P^''.,= (F'^'^Lj. + F^"^Lg).,-0. 



(63) 



We shall consider a spacetime metric representing a dynamic spherically symmetric (3 + 1)— dimensional wormhole 
given by ^. The non-zero Einstein tensor components arc the same that in Eq. ([7]). 

The electromagnetic tensor, compatible with the symmetries of the geometry, is given by ([5|). The F invariant is 
given by Eq. ^ while the invariant G is 
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The components of the stress-energy tensor, Eq. (I62p . in the orthonormal frame take the following form 
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^ft = r?i = (withi^j). (64) 

As in [l3|, Tjj = —Tfr, so the solution to Gj£ ~ —Gff can be solved separating variables like in (fTTj) : 

fe(r) = r [1 - a^ (^^2 _ ^2-)j ^ (■gg^ 

where a is a constant, and Ci and C2 are constants of integration. 
The electromagnetic field Eqs. (|551) take the form 



1-WM-^^ =0, (67) 



Eq. (pT)) can be solved for ELp 
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£;Lf = ^^ + ^1 (r) . (70) 
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Now, Eq. ([S5|) . when we take into account Eq. ([70]) give us the expression for Fi(r) 

Fi{r) = -g , C£=cte. 



Then ([701) takes the form 
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From this last relation we verify that ELp is singular at the throat. Integrating Eq. ([69)) . we obtain B = qmSnv{9). 
Of course when Lg = 0, we recover the expression obtained in [ij] for ELp. 

Therefore, when there is dependence on both electromagnetic invariants, F and G, we recover the solution in [17| . 
except for an additional term kG, 

L = L{F) + kG 

This term breaks the duality rotation symmetry of the electromagnetic field, present for instance in nonlinear 
electrodynamics of the Born-Infeld type. 

We conclude that geometry (4) does not allow nonlinear electromagnetic matter related to a lagrangian L{F,G). 
There is open still the question if assuming a less symmetrical spacetime it would admit such kind of matter. 

V. FINAL REMARKS 

In this article we have analyzed the solution representing an evolving wormhole coupled to nonlinear electrodynamics 
given by Arellano and Lobo in [17]. We presented the integration of the geodesic equations for light rays and the 
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corresponding graphics that show the influence of the nonhnear field; the traversability conditions as well as the weak 
field limit were analyzed. 

In the conclusions about traversability, we obtained that the tidal gravitational forces are constant in the radial 
movement, and they can be adjusted to allow a safe human passage. 

For the Maxwellian limit we have found that when r — > cxd the spacetime is asymptotically flat except for a remanent 
energy, a possible interpretation can be the presence of a magnetic perfect fluid that can be obtained from the stress 
energy tensor taking that limit. We must stress that the wormhole solution is only valid in a certain region that, even 
if it can be extended to be large, leaves outside the limit r —^ oo. However, the allowed trajectories for photons show 
that the NLED effective geometry is less restrictive than the one for the wormhole. 

When the dependence on both invariants for the lagrangian is considered, it turns out that the electric-magnetic 
duality symmetry is lost. 

Acknowledgments 

A. V. B. A. acknowledges a Ph. D. grant by Conacyt-Mexico. N. B. acknowledges partial support by Conacyt- 
Mexico, project 49182-F. R. G-S. was partially financed by Conacyt-Mexico, project 49924-J. 



[1] Einstein, A. and Rosen, N., "The particle problem in the General Theory of Relativity", Phys. Rev. 48, 73 (1935). 

[2] Wheeler, J. A., "Geons", Phys. Rev. 97, 511 (1955). 

[3] Morris, M. S., and Thorne, K. S. "Wormholes in spacetime and their use for interstellar travel: A tool for teaching General 

Relativity", Am. J. Phys. 56, 395 (1988). 
[4] Visser, M., Lorentzian Wormholes (American Institute of Physics, New York, 1995). 
[5] Lobo, F. N. Ph. D. thesis, Centro de Astronomia e Astrofisolutions in General Relativity: Traversable wormholes and 'warp 

drive' spacetimes", invited chapter to appear as edited collection "Classical and Quantum Gravity Research Progress", 

Nove Science Publishers [arXiv:gr-qc/0710.4474vl]. 
[6] Roman, T. A. "Inflating Lorentzian wormholes," Phys. Rev. D47, 1370 (1993) arXiv: gr-qc/9211012] . 
[7] Kar, S. "Evolving wormholes and the energy conditions," Phys. Rev. D49, 862 (1994). 
[8] Kar, S. and Sahdev, D. "Evolving Lorentzian wormholes," Phys. Rev. D53, 722 (1996) [arXiv:gr-qc/9506094| . 



[9] Anchordoqui, L. A., Torres, D. F. and Trobo M. L., "Evolving wormhole geometries," Phys. Rev. D57, 829 (1998) 

rarXiv:gr-qc/9710026; . 
[10] Kim, S. W. "Cosmological model with a traversable wormhole," Phys. Rev. D53, 6889 (1996). 
[11] Seiberg, N. and Witten, E., "String theory and noncommutative geometry," JHEP 9909, 032 (1999) 

arXiv:hep-th/9908142 . 
[12] Novello, M., Bergliaffa, S. E. P. and Salim, J. "Nonlinear electrodynamics and the acceleration of the universe," Phys. Rev. 

D69 127301 (2004) l arXiv:astro-ph/0312093| . 

Moniz, P. V. "Quintessence and Born-Infeld cosmology," Phys. Rev. D66, 103501 (2002). 

Garcia-Salcedo, R. and Breton, N. "Born-Infeld Cosmologies," Int. J. Mod. Phys. A15, 4341-4354 (2000) 
arXiv:gr-qc/0004017 . 

Garcia-Salcedo, R. and Breton, N. "Nonlinear electrodynamics in Bianchi spacetimes," Class. Quantum Grav. 20, 5425- 

5437 (2003) (arXiv:h ep-th/02 12T30] . 

Garcia-Salcedo, R. and Breton, N. "Singularity-free Bianchi spaces with nonlinear electrodynamics," Class. Quantum Crav. 

22, 4783-4802 (2005) ar Xiv:gr-qc /0410142 . 

Dyadichcv, V. V., Gal'tsov, D. V., Zorin, A. G. and Zotov, M. Yu. "Non-Abelian Born-Infeld cosmology," Phys. Rev. 

D65, 084007 (2002) arXiv:hep-th/0111099 . 

Vollick, D. N. "Anisotropic Born-Infeld Cosmologies," Cen. Relativ. Crav. 35, 1511-1516 (2003) arXiv:hep-th/0102187j. 
[13] Cataldo, M. and Garcia, A. "Non-singular charged black hole solutions for non-linear source," Cen. Relativ. Crav. 31, 

629-633 (1999) arXiv:gr-gc/9911084 . 

Burinskii, A. and Hildebrandt, S. R. "New type of regular black holes and particlelike solutions from nonlinear electrody- 
namics," Phys. Rev. D65, 104017 (2002) arXiv:hep-th/0202066; . 

Breton, N. "Stability of nonlinear magnetic black holes," Phys. Rev. D72, 044015 (2005) arXiv:hep-th/0502217\ 
[14] Bronnikov, K. A. "Regular magnetic black holes and monopoles from nonlinear electrodynamics," Phys. Rev. D63, 044005 

(2001) [arXiv:gr-qc/0006014| . 
[15] Baldovin, F., Novello, M., Perez Bergliaffa, S. E. and Salim, J.M. "A nongravitational wormhole". Class. Quantum Crav. 

17, 3265-3276 (2000) [arXiv:g r-qc/0003075 . 

Bronnikov, K. A. and Grinyok, S. V. "Electrically charged and neutral wormhole instability in scalar-tensor gravity", Crav. 

Cosmol. 11, 75-81 (2006) arXiv:gr-qc/0509062 . 



13 

[16] Arellano, A. V. B. and Lobo, F. S. N. "Non-existence of static, spherically symmetric and stationary, axisym- 

metric traversable w ormholes coupled to nonlinear electrodynamics", Class. Quantum Grav. 23, 7229-7244 (2006). 

|arXiv:gr-qc/0604095] . 
[17] Arellano, A. V. B. and Lobo, F. S. N. "Evolving wormhole geometries within nonlinear electrodynamics". Class. Quantum 

Crav. 23, 5811-5824 (2006). 
[18] Plebariski, J. F. , "Lectures on non-linear electrodynamics," monograph of the Niels Bohr Institute Nordita, Copenhagen 

(1968). 

Gutierrez, S. A., Dudley, A. L. and Plebahski, J. F. "Signals and discontinuities in general relativistic nonlinear electro- 
dynamics," J. Math. Phys. 22, 2835 (1981). 
[19] Lemos, J. P. S. , Lobo, F. S. N. and de Oliveira, S. Q. "Morris- Thorne wormholes with a cosmological constant," Phys. 

Rev. D 68, 064004 (2003) ■arXiv:gr-qc/0302049 . 
[20] The graphics were done using Dynamics Solver by J. M. Aguirregabiria, freeware program to analyze dynamical systems, 

available at <http://tp.lc.ehu.es/jma.html>. 
[21] Misner, C. W., Thorne, K. S. and Wheeler, J. A., "Gravitation," (Freeman, San Francisco, 1973), Sec. 32.6. 
[22] Bardeen J., presented at GR5, Tiflis, USSR, published in conference proceedings, 1968. Ayon-Beato E. and Garcia A., "The 

Bardeen model as a nonlinear magnetic monopole," Phys. Lett. B493, 149-152 (2000) arXiv:gr-qc/0009077 . 
[23] Dymnikova L, "Regular electrically charged vacuum structures with de Sitter center in Nonlinear Electrodynamics coupled 

to General Relativity," Class. Quantum Grav. 21, 4417-29 (2004) arXiv:gr-qc/0407072 . 

Burinskii, A. and Hildebrandt, S. R. S. R. "New type of regular black holes and particlelike solutions from nonlinear 



electrodynamics," Phys. Rev. D 65, 104017 (2002) arXiv:hep-th/0202066 



[24] Novello, M., Perez Bergliaffa, S. E. and Salim, J. M., "Singularities in general relativity coupled to nonlinear electrody- 
namics". Class. Quantum Grav. 17, 3821 (2000). 



